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Abstract 

In the paper regularity of solutions to stochastic Volterra equations in a separable 
Hilbert space is studied. Sufficient conditions for the temporal and spatial regularity 
of stochastic convolutions corresponding to the equations under consideration are 
provided. The results obtained generalize some well-known regularity results for 
solutions to stochastic differential equations. The paper is a continuation of previous 
author's papers concerning stochastic Volterra equations. 



1 Introduction 



Assume that (fi, J-", {J^t)t>o,P) is a stochastic basis. We study stochastic linear Volterra 
equations in a separable Hilbert space {H, | ■ of the form 

X{t) =Xo+ [ a{t- t) AX{T)dT + Wit) , t>0. (1) 
Jo 

In ([T]), Xq is an if- valued, J-Q-measurable random variable, a scalar kernel function a G 
Lj'Qj,(M_i_; M) and A is a closed unbounded linear operator in H with a dense domain 
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D{A). In the domain D{A) we introduce the graph norm | ■ \d{A) of A, i.e. |/i|D(A) : = 
+ lAhl'jjY^'^. In the paper we study the equation ([T]) which is driven by a Q- Wiener 
process W with TrQ < +00. 

Let us emphasize that the equation ([1]) includes a big class of equations, for instance 
heat and wave equations and integrodifferential equations. Moreover, it is an abstract 
version of several technical problems. 

Regularity of solutions to equations is very important and plays a prominent role in the 
study of evolution equations, particularly stochastic ones, see e.g. 
[ll[T5]. In the paper we expose sufficient conditions for time and space regularity of 
stochastic convolution, which is a solution to the stochastic equation ([1]). In our study we 
use the so-called resolvent approach to the equation ([1]), which is a natural extension of 
the semigroup approach used to stochastic differential equations. We note that the results 
obtained in the paper are an extension of the regularity results for solutions to stochastic 
differential equations given, e.g., in 

Let us note that our case considered in this paper is more difficult than the semigroup 
case connected with stochastic differential equations. The convolution appearing in ([1]) 
causes several serious problems and difficulties; some of them have been discussed in [8]. 
The most big problem comes from the fact that the solution operator (or the resolvent 
operator) corresponding to ([T]) and recalled below, does not create, in general, any semi- 
group. In consequence, powerful semigroup tools can not be used to the study of the 
equation ([1]). Particularly, the factorization method due to G. Da Prato, S. Kwapieh and 
J. Zabczyk, provided in [3], is not available in the resolvent case. However, some regu- 
larity results for solutions to ([T]) are known, e.g., those due to Ph. Clement and G. Da 
Prato [2J. 

Our approach to regularity of solution to ([I]) is different. We rewrite the stochastic 
convolution associated with ([1]), which is a mild and/or weak solution to ([1]). Next, we 
use the new form (see formula (|T2l) below) of the stochastic convolution and join it with 
appropriate Cauchy problem (see formulas f|T3l) and f|T5|) below). Then we use regularity 
results of some Cauchy problems considered in monographs [T2HT1], adapting them for 
our purposes. 

The paper is organized as follows. Section 2 contains an appropriate deterministic 
background while section 3 provides the main definitions and concepts used in the paper. 
In section 4 we give sufficient conditions for stochastic convolution to be regular solution 
to the equation ([T]). The final section 5 consists of proofs of the main results. 
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2 Deterministic Volterra equation 



The following equation 




t > 0, 



(2) 







is a deterministic counterpart of ([T]) in the space H. In the equation ([2]), the operator 
A and the kernel function a are the same as previously and / is a continuous if-valued 
function. 

In the paper we assume that the equation ([2]) is well-posed, that is, the family Sit), t > 
0, of the resolvent operators corresponding to the Volterra equation ([2]) exists and is 
defined as follows. 

Definition 1 (See, e.g., fTSf.) A family {S{t))t>o of bounded linear operators in H is 
called resolvent for the equation ^ if the following conditions are satisfied: 

1. S{t) is strongly continuous on and S{0) = I; 

2. S{t) commutes with the operator A, that is, S{t){D{A)) C D{A) and AS{t)x = 
S{t)Ax for all x G D{A) and t > 0; 

3. the following resolvent equation 



holds for all x e D{A), t>0. 

Because the class of locally integrable scalar kernel functions is too wide for obtaining 
good enough results concerning regularity of solutions to the equation ([1]), we shall con- 
sider the case when the kernel function is completely positive understood in the following 
sense. 

Definition 2 We say that function a G L^([0,T];M) is completely positive on [0,T], 
T < +00, if for any fi > 0, the solutions of the equations 




(3) 




and 





are nonnegative on [0,T]. 



There are several examples of completely positive kernels, e.g.: 
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1. a(t) = r-Vr(a), aG(0, 1], t > 0; 

2. a(t) = e-\ t > 0. 

The class of completely positive kernels appears in the theory of viscoelasticity. Several 
properties and examples of such kernels can be found in p[6;, Section 4.2]. 

In the paper we shall use the following theorem providing convergence of the resolvents 
for the equation ([2]). 

Theorem 1 ( Theorem 4]) Let A be the generator of Co-semigroup in H and suppose 
the kernel function a is completely positive. Then the pair {A, a) admits an exponentially 
bounded resolvent S{t), t > 0. Moreover, there exist bounded operators An such that 
{An, a) admit resolvent families Sn{t), t>0, satisfying ||S'„(t)|| < Me""'* (M > 1, wo > 
0) for all t > and 

Sn(t)x — )■ S(t)x as n — )■ +oo, (5) 

for all X & H, t > 0. Additionally, the convergence is uniform in t on every compact 
subset o/M_|_. 

In the above theorem and in the sequel, we use the Yosida approximation of the 
operator A and we denote it by 

An:=nAR{n,A) = n^R{n,A) -nl, neN. (6) 

The analogous results like Theorem [1] hold in other cases; for more details concerning 
the approximation results we refer to [7l[T0|[TT]. 

3 Probabilistic background 

Let Q G L{H) be a linear bounded symmetric nonnegative operator acting in the space 
H. We recall that the Wiener process W with the covariance operator Q is assumed to 
be a genuine one, that is, TrQ < +oo. 

In the paper we shall consider the following types of solutions to the equation ([T]). 

Let A* denote the adjoint of A with a dense domain D{A*) C H equipped with the 
graph norm | ■ \d{a*)- 

Definition 3 An H-valued predictable process X{t), t G [0,T], is said to be a weak 
solution to (Qy, ^f P{Jq \a{t - T)X{r)\HdT < +oo) = 1 and if for all ^ G D{A*) and all 
t G [0, T] the following equation holds 

{X{t),OH = {Xo, Oh + ( r a{t - r)X(r) dr, A*Oh + {W{t),OH, P-a.s. 

^0 
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As we have already written we assume that ([2]) admits a resolvent family S{t), t > 0. 
Then we are able to use the resolvent approach to both equations (E]) and ([I]). So, we can 
introduce the following definition. 

Definition 4 An H-valued predictable process X{t), t G [0,T], is said to be a mild 
solution to the stochastic Volterra equation if for arbitrary t G [0,T], 

X{t) = S{t)Xo+ [ S{t-T)dW{T), P-a.s. (7) 
Jo 

where S(t), t>0, is the resolvent for the equation 

We introduce the stochastic convolution, which is the main part of the mild solution 

W^{t):= [ S{t-T)dW{T), t>0. (8) 
Jo 

Let us recall some results concerning the convolution W^{t), t > 0, proved in [6]. 

Proposition 1 (See Propositions 3 and 4]-) Assume that ^ admits the resolvent 
operators S(t), t > 0. Then the process W^{t), t > 0, given by ^ has a predictable 
version. Moreover, the process W^{t), t > 0, has square integrable trajectories. 

In some cases weak solution of the equation ([1]) coincides with mild solution of ([1]); 
for more information we can refer to |6[[^. 

Proposition 2 (See JBi Proposition 5].) Let a G BV{M.^; M) and suppose that ^ admits 
a resolvent family S G C"^(0, oo; L(iJ)). Let X he a predictable process with integrable 
trajectories. Assume that X has a version such that P{X{t) G D{A)) = 1 for almost all 
t G [0,r]. If for anyte [0,T] and ^ E D{A*), 

{X{t),OH={Xo,OH + j\a{t-T)X{T),A*OHdT (9) 

+ {W{t),OH, P-a.s., 

then ^ 

X{t) = S{t)Xo+ [ S{t-T)dW{T), tE[0,T]. (10) 
Jo 

Proposition 3 (See JE, Proposition 6].) Assume that A is a closed linear unbounded op- 
erator with the dense domain D{A), a E Lj'Q^(R_|_; M) and S{t), 
i > 0, are resolvent operators for the equation Then the stochastic convolution 
fulfills the equation ^ with Xq = 0. 
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Corollary 1 (Compare ^ Corollary 3.10].) By Proposition\^ if A is a linear hounded 
operator, we have 

W^{t)= [ a{t-T)AW^{T)dT + W{t), te[0,T]. (11) 
Jo 

Let us note that the above formula (fTTj) indicates continuity of trajectories of the 
convolution W''^{t), t e [0,T], in the case when A is a bounded operator. In our paper 
we shall show continuity of trajectories in a more general case when A is an unbounded 
operator. Additionally, we will prove some regularity results for the convolution ([8]). 



4 Time and space regularity 



In this section we study the temporal and spatial regularity of the stochastic convolution 
W'^, defined by (ED, that is 

ft 

T < +00. 



w^{t)= [ s{t-T)dW{T), te[o,r], 

^0 



Theorem 2 Assume that the operator A in the equation ^ is the generator of a Co- 
semigroup T(t),t G [0,T]. Let the kernel function a(t),t G [0,T], be completely positive 
and, additionally, a G Lii„,([0, T]; R). Let S{t), W{t) and W^{t),t G [0,T], he Uke above. 
Then the following formula holds 

W^{t) = cA[ T{t-T) [ a{T -a)W^a)da + cW{T) dT + W{t), (12) 
Jo Uo J 

t G [0,T], where c = a(0) is a constant. 

Examples of functions, fulfilling assumptions of Theorem [2l are: 

1. a(t) = e"*, t G [0,r]; 

2. a(t) = r-Vr(a), aG(0, 1], te[e,T], e > 0. 



Comment One can see that the above formula (IT2l) is more complicated than the 
analogous one for the semigroup case, see, e.g. jH Chapter 5]. In fact, if a{t) = 1, the 
formula (|T2ll reduces to 

W^{t)=A[ T{t-s)W{s)ds + W{t), te[0,T]. 
Jo 



For our convenience we will assume in the sequel that c = a(0) = 1. 
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Corollary 2 Suppose that the assumptions of TheoremlE hold and a(0) = 1. Define the 
process 



where 



Then 



Y{t):= / T{t-s) W{s) + W{s) ds, t E [0,T] 



d{s - a) (a) da, sg[0,T]. 



(13) 



W{s) :-- 



W'^it) = AY{t) + W{t), te[0,T]. 



Additionally, Y belongs to C^{[0,T]; D{A)), F — a.s., and 
dY{t) 



dt 



AY{t)+ W{t) + W{t) 



tG [0,T]. 



(14) 



(15) 



From now we assume that the space if is a complex separable Hilbert space and the 
operator A is the generator of a strongly continuous analytic semigroup. 

Because both processes Y and W have continuous trajectories, directly from the for- 
mula (fT4l) we can deduce the following result. 

Theorem 3 Let us take the same assumptions like in Theorem and additionally let 
the semigroup he analytic. Then the stochastic convolution W^{t), t > 0, has continuous 
trajectories. 

Comment Let us note that Theorem [2], Corollary [2] and Theorem [3] can be formulated 
for other cases when the convergence ([5]) of the resolvents for the equation ([2]) holds. For 
instance in the cases described in the papers [71191 flO]. 



Unfortunately, we are not able to obtain Holderianity of trajectories of the stochastic 
convolution ([8]). In contrary to the semigroup case, the formula f|T2|) for the stochastic 
convolution is more complicated and contains the term W (see formula f|T3|) ). In conse- 
quence, the nonhomogeneous Cauchy problem ( fTSi) contains two functions W and W . The 
first one has Holder-continuous trajectories, but the second one has not such trajectories. 

There are available some temporal regularity results for stochastic Volterra equations 
obtained under different assumptions, see e.g., p] and [9]. 

In order to study next properties of the stochastic convolution ([8]) we shall introduce 
some scales of subspaces of H. We will use some results for analytic semigroups generated 
by sectorial operators and the Cauchy problems associated with such semigroups. 
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Definition 5 (See, e.g., JT^\T^ .) A linear operator A : D{A) C H ^ H is called 
sectorial if there are constants u eM., 9 E {71/2, n), M > such that 

(1) p{A)DSe,^:={\eC:\j^u, | arg(A - < ^} , 

(2) ||/2(A,A)|U(H)<-^ ^^^^ 

I A — uj\ 

Assume that the operator A is sectorial, then we can consider an analytic semigroup 
T{t), t > 0, defined in H through the Dunford integral. For details on sectorial operators 
and analytic semigroups we refer to the great monographs [HE] and [T2HI1]. If the 
semigroup T(t), t > 0, is assumed to be of negative type, we can define the fractional 
powers of operators as follows. For any 7 G (0, 1), we set 

i-A)-''x = —[ i-\)-miX,A)xd\, xeH, 
2vr2 Jr. , 

where A is a complex number and -R(A, A) := {XI—A)~^ is the resolvent operator of A. The 
curve r^e is defined as follows. For 6 G (7r/2,7r), 

r > 0, TrM = —Tig — r,^g + r^g, where Tlg,T^g are the half lines parametrized re- 
spectively hj z = C,e^^, z = C,e~^^, > r, and T^g is the arc of circle prametrized by 
z = re'"^, -6 <ri <e. 

By {—A)'^ we shall denote the inverse of the operator {—A)~'^ and by D{{—A)'^) its 
domain. In the literature the operators {—Ay are called the fractional powers of the 
operator —A. They have very interesting and useful properties, see e.g. [T2HII]. Among 
others, there exist contants M^, M^^, /c = 0, 1, 7 G (0, 1), such that 

||/i^r(t)|| < Mfc^-^ t>o, (17) 

and 

1 1 i-AyA'' T{t)\\< Mk^ t-^~\ t > 0. (18) 

In the remaining part of the paper we shall use the following properties of powers of 
operators. 

Theorem 4 (IM Theorem 4.6]) Let a,(3 eC such that Re /3 < Re a. Then D{A"') C 
D{A^), and for every x G D{A°') A^x = A^^"-A"-x. Moreover, for each x G D{A^), 
A^x G /^(A^"") and A^'^A^x = Conversly, if x e D{A^) and A^x G D^A'^'P), 

then X G D{A°') and again A^'^'^A'^x = A'^x. 

Now, we are able to formulate spatial regularity for the convolution ([8]). 
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Theorem 5 For all t G [0,T] and 7 G (0, 1), the convolution W'^{t) belongs to the space 
D{{—Ay), a.s., provided the process W belongs to the domain D{{—A)'^). Additionally, 
{—AY'W^ is a Gaussian process and for all 

7 G (0, 1) trajectories of the process {—A)'^W^ are continuous. 

Following A. Lunardi [l2l[T3] we introduce the interpolation spaces. 

Definition 6 For 7 G (0, 1), I <p < +00, we denote Da{j,p) '■= {H, D{A))^^p, where 
the pair {H, D{A))^^p is the interpolation spaces according to e.g. Definition 1.2 in jl3^ . 

Definition 7 For 7 G (0, 1), 1 < p < +00, we define Da{i + := {x G D{A^) : 
A^x G D{A)^^p\, with the norm IxlD^i-y+k^p) '■= \x\h + l^^xloAi-y,?)- 

Let us note that -0^(7 + k,p) is the domain of the part of A^ in Da(7,p)- 
Now, we are able to formulate further regularity results for the stochastic convolution 
defined by ([8]) in terms of interpolation spaces. 

The formula f|T5l) enables us to study regularity of the process Y using properties of 
deterministic Cauchy problem. 

By [ini Proposition 6.6] and from square integrability on [0,T] of trajectories of the 
processes W and W we have the following result. 

Theorem 6 For every 7 G (0,1), trajectories of the process Y belong to the space 
L2((0, T); L)^(7, 2))niy^'2((0, T);H). Consequently, it belongs to the space W^-^''^{{0, T); 2)) 
for every 9 < •y < 1 and to C([0, T]; D^(7, 2)) for every 7 G (0, |). Moreover, there is a 
constant M independent of the processes W and W , such that 

\Y\w'y''^{(0,T);H) + |^|l2({0,T):Da(7,2)) < M \W + | L2((0,T);/f) • 

The next result is formulated without using interpolation of powers of operators. We 
shall use so-called Rademacher bounded operator. 

Definition 8 A subset of L{H) is W-bounded (Rademacher bounded) if there is C > 
such that for all n & N, Li, . . . , L„ G L, a;i, . . . , x„ G if we have 



n 

rkTkXk 


< 


n 

rkXk 




k=l 




k=l 


L2({0,1:H)) 



where rk{t) := sign{sm2^Tit) are Rademacher functions in (0, 1). 
Using Theorem 6.11] we can formulate the following theorem. 
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Theorem 7 Assume that the operator A 
group generated by A is strongly continuous, 
alent: 



D{A) G H ^ H is sectorial and the semi- 
Then the following two statements are equiv- 




process 



Y 



belong 



to 



the 



space 



2. there are u 



e 



e 



> 



such that the family of operators 



{(A — u)R{\, A) : \ E Sg^^} is IZ-bounded. 

Let us note, see e.g. remark on page 171 in [13], that in Hilbert space every sectorial 
operator satisfies the condition (2) in Theorem [71 So, we can formulate the following 
conclusion. 

Corollary 3 Assume that A : D{A) d H ^ H is a sectorial operator and the semi- 
group generated by A is strongly continuous. Then trajectories of the process Y belong to 
W''\{0,Ty,H)nL\iO,Ty,D{A)). 

5 Proofs of the main results 

Proof of Theorem [2l 

Because the formula (fTT!) holds for any bounded operator, then it holds for the Yosida 
approximation An of the operator A, that is 



^0 

In the formula ffH?]) , Snit), t > 0, is the resolvent genereted by the pair {A^, a{t)), t >0, 
(see Theorem [T]) and 




(19) 




(20) 



Let us denote 




(21) 



Then, from the Leibniz rule 




(22) 
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From (1191) and ( 1211) . we can write 

W'"{t)=A^Z„.{t) + W{t), tG[0,T] 
From ([22]), if a(0) ^ 0, we have 

1 



a(0) 



Z'^{t)- / a(t-r)Vr^"(r)c/r 



So, we obtain 



Z'^{t)- [ a{t-r)W'-{r)dr = a{0)[AnZ4t) + W{t)], t e [0,T]. 
Jo 



And next 



Z'^{t) = a{0)AnZn{t)+ [ a{t-T)W^"{T)dT + a{0)W{t), tG[0,T]. 

Jo 



For simplicity, we set 



W^{t):= / d(t-r)Vr^"(r)dr, te[0,T]. 



So, we have 

Z'^{t) = cAn Zn{t) + [Wf (t) + cW{t) \ , where c = a(0). 
Prom the above formula 

Znit) = [ e^(*-^)^" \w^{t) +cW{t)] dr, t G [0,T] , 

because ^„(0) = 0. 

From the formula (1T9|) . 

= A„Z„(t) + H^(t), tG[0,T], or 
= AJ„Z„(t) + t G [o,r], 

where J„ := ?T,i?(?T,, A). 
So, we have 

W^^it) = AJn [ e"(*-")^" \w^{r) + cW{t)\ dr + W{t), t G [0, T]. 
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Basing on Theorem [T], properties of Yosida approximation An of the operator A, and the 
Lebesgue dominated convergence theorem, we have: 

hm JnX = X, for any x & H; 

n—^oo 

hm AnX = Ax, for any x G D{A); 

n— >oo 

hm e^^" X = T(t) x, for any x G H; 



and hm sup E W^^ (t) - (t) = 0. 

"^-^•^ te[o,T] 

Because the operator A is closed, we can conclude that the integral 
T{t-T) W^{t) + cW{t) dr belongs to the domain £)( A). 
Hence, passing to the limit with n — )■ +00, we obtain 

W^{t) = cA [ T{t-T) \w^{t) + cW{t)] dT + W{t), te[0,T], 
Jo ^ 



where 

W'^ir) = [ d(r - a)W^{a)da , r G [0, T] . 
Jo 

In the case a(0) = 0, from the formula ( 1T9|) . passing to the limit we would have only 

W^{t)= [ a{t-T)AW^{T)dT + W{t), te[0,T], 
Jo 

that is, the formula like (fTTj) . ■ 
Proof of Corollary O 

From the definition (formula ( 1T3|) ) of the process Y and properties of convolution, Y G 
C^{[0,T]; D{A)), P — a.s. Next, from the Leibniz rule and property of semigroup we 
obtain 



dY{t) 
dt 



* dT{t - s) 



W{s) + W{s) ds + T{0) W{t) + W{t) 



A / T{t-s) W{s) + W{s 



ds + 



W{t) + W{t) 



AY{t)+ W{t) + W{t) 



tG [0,T]. 



Proof of Theorem [5l 

In our case we have no an appropriate representation formula for the resolvent like in the 
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semigroup case, see, e.g., [TH Chapter 2]. In consequence, we have no analogous estimates 
for resolvents like (IT71)- (|T8l) . In order to omit these difficulties, we have to use the formula 
f[T^ for the convolution W^. 

By (|T8|) . the integral 



Jo '- 
So, the process {—AyW^ is well defined. 

Gaussianity of the process comes from Gaussianity of the process . Continuity of 
trajectories of the process {—AyW^ results from continuity of trajectories of the process 
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